Influence of magnetic fields on cold collisions of polar molecules 
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We consider cold collisions of OH molecules in the 2 n 3 / 2 grounds state, under the influnce of a 
magnetic field. We find that modest fields of several thousand gauss can act to suppress inelastic 
collisions of weak-field seeking states by two orders of magnitude. We attribute this suppression to 
two factors: (i) An indirect coupling of the entrance and the exit channel, in contrast to the effect of 
an applied electric field; and (ii) the realtive motion of the entrance and exit scattering thresholds. 
In view of these results, magnetic trapping of OH may prove feasible. 

PACS numbers: 34.20.Cf,34.50.-s,05.30.Fk 



I. INTRODUCTION 

As the experimental reality of trapping ultra-cold polar 
molecules approaches, a clear understanding is needed of 
how the molecules interact in the trap environment. On 
the most straightforward level, collisions are essential for 
cooling the gas by either evaporative or sympathetic cool- 
ing methods. A high rate of elastic collisions is desirable, 
while a low rate of exothermic, state-changing collisions is 
essential if the cold gas is to survive at all. Furthermore, 
a clear understanding of 2-body interactions allows one 
to construct a realistic model of the many body physics 
in this dilute system [lj. 

One promising strategy for trapping ultracold 
molecules might be to follow up on successes in trapping 
of cold atoms, and to construct electrostatic 0,0 or mag- 
neto static jj| traps that can hold molecules in a weak- 
field-seeking state. Cold collisions of polar molecules in 
this environment have been analyzed in the past, find- 
ing that inelastic collision rates were unacceptably high 
in the presence of the electric field, limiting the possi- 
bilities for stable trapping 0. Ref. Q found that the 
large inelastic rates were due to the strong dipole-dipole 
interaction coupling between the molecules. One impor- 
tant feature of the dipole-dipole interaction is its com- 
paratively long range. Even without knowing the details 
of the short-ranged molecule- molecule interactions, the 
dipole forces alone were sufficient to change the inter- 
nal molecular states. Indeed, a significant finding was 
that for weak-field seekers, the molecules are prevented 
from approaching close to one another due to a set of 
long-range avoided crossings. Therefore, a reasonably ac- 
curate description of molecular scattering may be made 
using the dipolar forces alone 0. 

A complimentary set of theoretical analyses have con- 
sidered the problem of collisional stability of paramag- 
netic molecules in a magnetostatic trap. For example, 
the weak-field-seeking states of molecules are expected 
to survive collisions with He buffer gas atoms quite well 
0,13- Collisions of molecules with each other are also ex- 
pected to preserve their spin orientation fairly well, and 
hence remain trapped [9(- However, this effect is miti- 
gated in the presence of a magnetic field 0, |H| . 

So far, no one appears to have considered the influ- 



ence of magnetic fields on cold molecule-molecule colli- 
sions where both species have electric dipole moments. 
In this paper we approach this subject, by considering 
cold OH( 2 n 3 /2)-OH( 2 n 3 /2) collisions in a magnetic field. 
To the extent that the applied electric field is zero, one 
might expect that dipole forces average to zero and thus 
do not contribute to de-stabilizing the spin orientation. 
It turns out that this is not quite correct, and that dipole- 
dipole forces still dominate long-range scattering. How- 
ever, applying a suitably strong magnetic field turns out 
to mitigate this effect significantly. Interestingly, even 
in this case the residual second-order dipole interactions 
are sufficiently strong to restrict scattering to large inter- 
molecular separation 

The main result of the paper is summarized in Figure 
^ which contrasts the influence of electric and magnetic 
fields. Figure [Ja) plots the elastic (solid curves) and 
inelastic (dashed curves) collision rate constants for OH 
molecules in their \FMpe) = |22— ), weak-field-seeking 
hyperfine state (for details on quantum numbes, see be- 
low). As the electric field is increased, the elastic rate 
constant grows to alarmingly large values, making the 
gas collisionally unstable, as was shown in Ref. Q. Fig- 
ure^b), which is new to this paper, shows the analogous 
rate constants in a magnetic field (in both cases the field 
is assumed to lie along the positive z axis of the labora- 
tory reference frame). In this case the magnetic field has 
the effect of suppressing collisions, all the way down to a 
rate constant of 2 x 1CT 13 cm 3 /sec at fields of B = 3000 
gauss. These results are moreover fairly robust against 
raising the temperature to the merely cold (not ultracold) 
temperatures, ~ lmK, attainable in buffer-gas loading 
or Stark slowing experiments. This is good news for ex- 
periments - it implies that cooling strategies that rely on 
collisions may be feasible, provided a suitably large bias 
magnetic field is applied. 

Our main goal here is to analyze the suppression of 
rates in a magnetic field. The organization is as follows. 
First we review the relevant molecular structure, and es- 
pecially the Stark and Zeeman effects, to illustrate their 
complementary natures. We then present an overview of 
the scattering model, including a review of the dipole- 
dipole interaction. Finally we present an analysis of the 
system in a magnetic field using a reduced channel model 
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FIG. 1: Thermally averaged rate constants for collisions of 
weak-field seeking sates of OH, as a fucntion of applied elec- 
tric field (a) and magnetic field (b). In both cases, solid lines 
denote elastic scattering rates, while dashesd lines denote in- 
elastic scattering rates. Two temperatures are considered. 
Applying an electric fields drives the inelastic collisions rates 
up while applied magnetic field drives inelastic collision rates 
down. 



that encapsulates the essential collision physics. Finally 
the model is qualitatively understood using the adiabatic 
representation. 



II. ELECTRIC VERSUS MAGNETIC FIELDS 
APPLIED TO MOLECULES 

Both the Stark and Zeeman effects in molecules have 
a similar form, since both arise as the scalar product of 
a dipole moment with an external field. Their influence 
on the molecule is quite different, however, since they act 
on fundamentally different degrees of freedom. The elec- 
tric field is concerned primarily with where the charges 
are in the molecule, whereas the magnetic field is con- 
cerned with where they are going. This is of paramount 
importance, since it implies that the electric field is a 
true vector (odd under the parity operation), whereas 
the magnetic field is a pseudovector (even under parity) 
[l2|. An electric field will therefore mix parity states of 
a molecule, while a magnetic field will mix states only 
of a given parity. This distinction is explored further in 
Ref. nere we wm om y focus on those aspects of 



immediate relevance to our project. 

The rest of this paper will, fundamentally, restate this 
fact in the context of scattering, and follow up the con- 
sequences that arise from it. To set the context of this 
discussion, and to fix our notation, we first consider the 
molecules in the absence of external fields. 



A. Molecular structure in zero external field 

The OH radical has a complicated ground state struc- 
ture which includes rotation, parity, nuclear spin, elec- 
tronic spin and orbital degrees of freedom. We assume 
that the vibrational degrees of freedom are frozen out at 
low temperatures, hence treat the molecules as rigid ro- 
tors. We further assume that perturbations due to far 
away rotational levels are weak. We do, however, include 
perturbatively the influence of the Q = 1/2 fine structure 
level, as described in Ref. 

The electronic ground state of OH is 2 n, with fi = 3/2 . 
OH is an almost pure Hund's case (a) molecule, meaning 
the electronic degrees of freedom are strongly coupled 
to the intermolecular axis. The electronic state of the 
molecule in the J basis is denoted by | JMjf2)|A£) where 
J is the total electronic angular momentum, Mj is its 
projection onto the lab axis and Q is J's projection onto 
the molecular axis. £ and A are the projection of the 
electron's spin and orbital angular momentum onto the 
molecular axis, and their sum equals £1 (= X + A). The 
electronic degrees of freedom, |AS), will be suppressed 
for notational simplicity because they are constant for 
all the collisional processes we consider. 

To describe the molecular wave functi on we assume 
a rigid rotor, (a,0,'y\JMjSl) = yf^D& jn (a, 0, 7 ), 

where(a, 0, 7) are the Euler angles defining the molecular 
axis and Djjf n is a Wigner D-function. It is necessary 
for a n state molecule to use the parity basis because OH 
has a A-doublet splitting which separates the two parity 
states (e//). The A-doublet arises from a coupling to a 
near by S state. It is the coupling of the X state to n 
state of the same parity which breaks the degeneracy of 
the two n parity states |l4j . 

In the parity basis the molecular wave function is writ- 
ten 



\JMjVt±) = 



\JMjil) + e\JMj - fi) 
71 



(1) 



Where e = +(— ) represents the e (/) state, and Ct = \ 
It should be noted that the sign of e is not the parity, 
rather parity is equal to e(— 1)( J_1 / 2 ) 14]. Thus for the 
ground state of OH where J = 3/2, parity is equal to 
— e. Throughout the paper we use ± denote the sign 
of e, not parity. In the parity basis there is no dipole 
moment, because this basis is a linear combination of 
electric dipole "up" and "down" . This fact has important 
implications for the dipole-dipole interaction. 

Including the hyperfme structure is important because 
the most dominant loss processes are those that change 
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the hyperfine quantum number of one or both of the scat- 
tering molecules. The hyperfine structure arises from 
interaction of the electronic spin with the nuclear spin 
(J) which must then be included in the molecular basis 
set. In the hyperfine basis the OH states are \FMp±), 
where F = J + 1 and Mp is its projection onto the lab 
axis. Here we suppress VL in the notation, as its value is 
understood. To construct basis functions with quantum 
number F we expand in Clebsh-Gordon coefficients: 

\FM F ±)= |JMjn±)|7Mj>(JMjJMr|FAfp>. 

MjMi 

(2) 

Relevant molecular energy scales for the scattering 
problem are the A-doublet splitting which is A ~ 
0.0797AT, the hyperfine splitting is A hf ~ 0.0038AT. 
OH also has an electric dipole moment is fi ~ 1.668D. 
Throughout this paper we use Kelvin (K) as the Energy 
unit, except in the instances of thermally averaged ob- 
servables. For reference, 0.64AT — 1cm -1 . 



B. Stark Effect in OH 

As noted above, the distinguishing feature of the Stark 
effect is that it mixes molecular states of opposite par- 
ity separated by the A doublet splitting. A consequence 
of this is that the Stark energies vary quadratically with 
electric field at low fields, and linearly only at higher 
fields. The field where this transition occurs is given 
roughly by equating the field's effect /2g • E to the A dou- 
blet splitting (here /2e is the molecule's electric dipole 
moment, and E is the field. In OH, this field is approxi- 
mately E ~ A/2fi E ~ lOOOlV/ cm).) 

The Stark Hamiltonian has the form 



Hs = -Jae ■ E 



(3) 



where we take the field to be in the z direction. In the 
basis in which f2 has a definite sign, the matrix elements 
are well known 1151: 



(JMjVt\H s \JMjVt) 



-/j, E EfLMj 
J(J + l) 



(4) 



In the Stark effect there is a degeneracy between states 
with the same sign of flMj, meaning ±Mj are degenerate 
in an electric field. We can recast the Stark Hamiltonian 
into the J-parity basis set from Eqn. (I). Doing so, we 
find 



(JMjne\H s \JMjCle') 



-fisEQMj (1 — ee 



J(J+ 1) 



(5) 



In this expression, the factor (1 — ee')/2 explicitly repre- 
sents the electric field coupling between states of opposite 
parity, since it vanishes for e = e' . 

Finally, using the definition of the F-parity basis in 
Eqn. P|l. we arrive at the working matrix elements of 



the Stark effect: 



(FM F e\H s \F'M F e') = - t i E E 



l + ee'(-l) 



,7+./'+2Q+l 



X(-1) 



J+J' +F+F' + I-M F -Q,+1 



[F,F',J,J'} 



F' 


1 


F 1 




F' 




Mp 





—Mp ) 




J 





.(6) 



In this notation [ji,j 2 ,...] = y/(2 3l + l)(2j 2 + 1)(...). 
Figure [21 shows the energy levels of OH in the presence 
of an electric field. Both parity states are shown, la- 
beled e and /. An essential point of Fig. [21 is that the e 
and / states repel as the electric field in increased. This 
means the all of the / (e) states increase (decrease) in 
energy as the field in increased, implying that states of 
the same parity stay close together in energy as the field 
is increased. This fact has a crucial effect on the inelastic 
scattering as we will show. 

The highest-energy state in Fig. [21 is the stretched 
state with quantum numbers \FMpe) — |22— ). It is this 
state whose cold collisions we are most interested in, be- 
cause i) it is weak-field seeking, and ii) its collisions at 
low temperature result almost entirely from long-range 
dipole-dipole interactions ,5j . Molecules in this state will 
suffer inelastic collisions to all of the other internal states 
shown. The rate constant shown in Figure ^ is the sum 
of all rate constants for all such processes. 



C. Zeeman Effect in OH 

When OH is in an external magnetic field the elec- 
tron's orbital motion and intrinsic magnetic dipole mo- 
ment both interact with the field. The interaction is de- 
scribed by the Zeeman Hamiltonian which is 



H z = -[is ■ B = ^ (L + g e S) ■ B. 



(7) 



Here /iq is the Bohr magneton and g e is the electron's g 
factor (<? e ~ 2.002). As above, we assume the field to be 
in the laboratory z direction. In the J basis, the Zeeman 
Hamiltonian takes the form |l6j 



(jMjn\H z \JMjn) 



(8) 



J(J + 1) 

This is quite similar to the equivalent expression for 
the Stark effect, except that the electron's g-factor plays 
a role. Interestingly, for a 2 H state the prefactor (A + 
g e E)0 is always greater than zero. We now recast the 
Zeeman interaction into the J-parity basis set QJ. This 
gives us 

VoB (A + q e t) VLMj 
(JMjne\H z \JMjne') = w V -5 te , (9) 



J(J + 1) 



for £1 = 3/2 states, and 



(JMjSle\H z \JMjne') = ~j(T^i) 5 ~> (10) 
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FIG. 3: Zeeman effect for the ground state of OH, in low (a) 
and high (b) fields. This plot is the same for both the e and 
/ states for zero electric field because the Zeeman interaction 
respects preserves parity. The [22—) state is indicated in gray. 



copies of the same energy level diagram, separated in en- 
ergy by the lambda doublet energy. For small magnetic 
fields the molecular g-factor is g^ag ^ {F 2 + J 2 — I 2 ), 
and is always positive for OH. This is in contrast to 
the low field magnetic moment of alkali atoms which is 
9ma<j % (F 2 — J 2 — I 2 ) (and where J, of course, refers to 
the sum of orbital and spin angular momenta). In Eqn. 
for Q = 1/2 the factor A + g e T, goes to A - g e T,. 



FIG. 2: Stark effect for the ground state of OH with the 
hyperfine structure accounted for. In zero field the f states 
and the e states are separated by the A-doublet energy. The 
Gray line indicates the state of interest for our analysis, the 
22—) state. An important feature of this interaction is that 
the opposite parity states repel and thus like parity state stay 
close together in energy. 



for = 1/2 states. Notice that for f2 = 1/2, the orbital 
and spin contributions to the molecular magnetic mo- 
ment nearly cancel, to within the deviation of ,g e /2 from 
one. For the f2 = 3/2 states of interest to us, however, 
the magnetic moment remains large. 

The key feature of the Zeeman matrix element © is 
that it is diagonal in e, in contrast to the Stark matrix 
element. This trait persists in the hyperfine basis as well, 
where the matrix elements are 



(FM F e\H z \F'M F e') = 

- /l + ee'(-D (J+J ' +2n) \ 
M0 5(A + 9eS) I 1 2 J J 

x ^yl+J'+F+F'+I-Mr-n+l ^ ^/^ ^ j/j 

J I J' 
~n o fi' 



(11) 



F' 


1 


F \ 




F' 




M F 





-M F ) 




J 





Figure |21 shows the Zeeman energies in the hyper- 
fine basis, for low (Fig. Efa)) and high (Fig. EJb)) 
fields. For OH in the 2 n 3 / 2 state, the parity factor 

/ ia- ee 'f_i')( J + J ' + 2fi ) \ 

I — — - — ^ I reduces simply to S ee i . Because the 

magnetic field respects parity, figure|21(b) amounts to two 



III. SCATTERING HAMILTONIAN 

A complete potential energy surface for the interac- 
tion of two OH molecules, including the relatively long- 
range part most relevant to cold collisions, is at present 
unavailable. Certain aspects of this surface have, how- 
ever, been discussed [13. For the time being, we will 
follow our previous approach of focusing exclusively on 
the dipole-dipole interaction. It appears that molecules 
in the highest weak-field-seeking states are mostly insen- 
sitive to short-range effects. 

The "raw" scattering channels have the form 
\FiM Fl ei)\F2M F2 e2)\lmi), which specifies the internal 
state of each molecule and the partial wave / describing 
the relative orbital angular momentum of the molecules. 
In a field, of course, the hyperfine and parity quantum 
numbers are no longer good. It is therefore essential to 
consider a set of scattering channels "dressed" by the 
appropriate field. This is achieved by diagonalizing the 
Stark or Zeeman Hamiltonian of each molecule, includ- 
ing the A doubling and hyperfine structure. The result- 
ing eigenvectors then comprise the molecular basis used 
to construct the scattering Hamiltonian. Field dressing 
is essential because otherwise non-physical couplings be- 
tween channels presist to infinite separation. The diago- 
nal contributions of the Stark and Zeeman Hamiltonian 
in the field dressed basis define to the scattering thresh- 
olds as R — > 00 . 

The channels involved in a given scattering process are 
further constrained by symmetries. Namely, scattering of 
identical bosons restricts the basis set to even values of I 
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only. In addition, the cylindrical symmetry enforced by 
the external field guarantees that the total projection of 
angular momentum on the field axis, M = Mp 1 + Mp 2 + 
mi, is a conserved quantity. 

The scattering wave function is expanded in these field- 
dressed channels, leading to a set of coupled-channel 
Schrodinger equations 



-ft 2 d 2 



2/2 



-1 



2m r dR 2 - 2m r R 2 



1 + V(R) + H FS ) $(R) = E$(R), 



(12) 

where ip is the multichannel wavefunction and m r is 
the reduced mass. The operators FLp$ denotes the fine 
strucutre, including the effect of an electric or magnetic 
field. (In this paper we do not yet include the simulta- 
neous effect of both fields.) 

In Eqn. (|12(l . the operator V_ represents the interac- 
tion between the molecules. We are most interested in 
the dominate dipole-dipole interaction whose general op- 
erator form is 



3(R • Ai)( R • Ag) - Ai - A2 
i? 3 



(13) 



where fa the electric dipole of molecule i, R is the inter- 
molecular separation, and R is the unit vector defining 
the intermolecular axis. This interaction is convenie ntly 
re- written in terms of tensorial operators as follows [18j : 



H, 



(Ml ®^lf- Q 



(14) 



Here C 2 (6,(f>) is a reduced spherical harmonic that acts 
only on the relative angular coordinate of the molecules, 
while (/Ji <g> H2)- q is the second rank tensor formed from 
two rank-1 operators, that act on the state of the 
jth mo lecule. These first rank operators are written as 
reduced spherical harmonics, Cg(atf3) where a and (3 are 
two of the Euler angles of the rigid rotator wavefunction. 
With this form of the dipole-dipole interaction, we can 
then evaluate the matrix element. 

In the hyperfine parity basis J2J the matrix elements 
are [f| 



(MmilH^llWm'i) 



[1,1 , Ji, J 1 , J 2 , J 2 , ri, f 1 , t 2 , t , 2 \ 



+2n 2 +i 



^3 L"ri"i5"l!-ai-ii!-i!-l)-i)-ZJ I 2 

X ^_l^(l+F l +F[+F-2+F^ + J 1 + J l 1 + J2 + .J' 2 +M l +M 2 -n' 1 -n'. 2 +Mi) 

Jl 1 J[ 



1 



1 



J 2 1 4 



M Fl - M F , M F2 - M F , Mi - M[ J \ - fit ^1 J \ -0 2 n 2 

1 F, F{ \f 1 F 2 

M Fl - M F{ -M Fl M F[ J \ M F2 - M F , ~M F2 M f , 

I' 2 I \ (V 2 I \ f F 1 F[ 1 \ f F 2 F' 2 1 

M v M t - Mi, -Mi) I J 1 J{ Jx I j 1 4 J 2 I 



(15) 



A central feature of this matrix element is the factor 
(1 - eiei) (1 - e 2 e' 2 ), using J { = J< = a t = 3/2. As a 
consequence of this factor, matrix elements diagonal in 
parity identically vanish in zero electric field. Instead, for 
example, two /-parity states only interact with one an- 
other via coupling to a channel consisting of two e-parity 
states. 

This dependence on parity is perhaps not unexpected, 
since the dipole-dipole force is of course transmitted by 
the dipole moment of the first molecule producing an 
electric field that acts on the second molecule. But in 
a state of good parity, the first molecule does not have 
a dipole moment until it is acted upon by the second 
molecule. Thus, both molecules must simultaneously mix 
states of opposite parity to interact. Notice that in the 
presence of an electric field, the dipoles are already par- 
tially polarized, and this restriction need not apply; the 



scattering channels are already directly coupled. This 
change is of decisive importance in elucidating the influ- 
ence of electric fields on collisions. In a magnetic field, by 
contrast, parity remains conserved and the interactions 
are intrinsically weaker as a result. 



IV. INELASTIC RATES OF OH-OH 
COLLISIONS IN EXTERNAL FIELDS 

We move now to the consequences of the interaction 
(|T5|l on scattering. Scattering calculations are done us- 
ing the log-derivative propagator method . To ensure 
convergence at all collision energies and applied fields, 
it was necessary to include partial waves up to Z = 6, 
and to carry the propagation out to an intermolecular 
distance of R = 10 4 (a.w.) before matching to long-range 
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wave functions. Cross sections and rate constants are 
computed in the standard way for anisotropic potentials 

We remind the reader that throughout we consider col- 
lisions of molecules initially in their \FMpe) — |22— ) 
states, which are weak-field seeking for both electric and 
magnetic fields. Thus for a scattering process incident on 
an s- partial wave, the incident channel will be written 
\i) = {FtMp^FiM^llmi) = |22,-)|22,-)|00). 

In the following, we will make frequent reference to 
"energy gap suppression" of collision rates. This notion 
arises from a perturbative view of inelastic collisions, in 
which case the transition probability amplitude is pro- 
portional to the overlap integral 

[ dRil>i(R)V if (R)il>f(R) (16) 

where ipij denote the incident and final channel radial 
wave functions, and Vif is the coupling matrix element 
between them. In our case, ipi will have a long de Broglie 
wavelength corresponding to its essentially zero collision 
energy. The de Broglie wavelength of ipf will instead 
grow smaller as the energy gap E- L — Ef between inci- 
dent and final thresholds grows. Thus the integral in 
(|16f) . and correspondingly the collision rates, will dimin- 
ish. For this reason, the collisions we consider tend to 
favor changing the hyperfine states of the molecules over 
changing the parity states, since the hyperfine splitting 
of OH is smaller than the A-doubling. 

A. Electric Field Case 

To calculate scattering in the presence of an electric 
field, we only need to include partial waves I = 0, 2 for 
numerical accuracy of K™ el for the field range that we 
consider, E < 1000(V/ cm), and at a collision energy of 
10~ 5 K. Here we are only interested in the trend and 
identification of the loss mechanism. To numerically con- 
verge the inelastic rates at higher field values, where the 
induced dipoles are large, naturally requires more partial 
waves. 

Figure0](a) shows the total (black) and partial (color) 
inelastic rate constant K™ el as a function of the electric 
field (compare Fig. 0(a)). Even in zero field, where the 
dipolar forces nominally average out, the rate constant is 
large, comparable to the elastic rate constant. This fact 
attests to the strength of dipolar forces in OH, even in 
second order. 

The green line in Fig. 0(a) represents losses to the 
dominant zero field loss channel |10— )|22— )|22). The 
blue curve in Fig. 01 (a) represents instead the domi- 
nant loss process at higher electric field values, in channel 
1 21 — ) 1 22 — )|21). Whereas the former rate remains rela- 
tively insensitive to field, the latter rises dramatically. 

This behavior arises from two competing tendencies in 
an electric field. The first is the increasing mixing of dif- 
ferent parity states as the field is turned on, leading to 




Electric Field (V/cm) 
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200 400 600 800 1000 
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FIG. 4: (a) Total (black) and selected partial (color) inelastic 
rates for OH-OH collisions as a function of electric field. The 
green curve is the dominant zero field inelastic loss process 
to channel |10— ) 1 22 — ) (green). In the presence of the field, a 
different channel, |21— ) 1 22 — ), becomes dominant (blue), (b) 
The thresholds for these exit channels, relative to the incident 
threshold. 

an increasing strength of the direct dipole-dipole coupling 
that affects both exit channels. This additional coupling 
would, in general, cause inelastic rates to rise. It is, how- 
ever, offset by the competing tendency for inelastic rates 
to become less likely when the change in relative kinetic 
energy of the collision partners is larger. Fig. 0(b) shows 
the threshold energies for the two exit channels in Fig. 
0(a), versus field, with zero representing the energy of 
the incident threshold. Here it is evident that loss to 
the channel |22— )|10— )|22) (green line) is accompanied 
by a large gain in kinetic energy, whereas loss to chan- 
nel 1 21 — ) 1 22 — ) |21) (blue line) gains comparatively little 
kinetic energy, and thus the later channel more strongly 
affected by the increased coupling generated by the field. 



B. Magnetic Field Case 

To gain insight into the suppression of the inelastic 
rates in a magnetic field (Fig. 0(a)), calculations were 
at a representative collision energy E = 10~ 5 K. To con- 
verge the calculations in high field (B> 1500 gauss) re- 
quired partial waves I = 0,2,4,6. We have only con- 
sidered collisions with incident partial wave / = 0, since 
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FIG. 5: (a) The total (black) and partial (color) inelastic 
rates for OH-OH collisions as a function of magnetic field. 
The colors are explained in the text, (b) The corresponding 
thresholds, referred to the incident channel's threshold (Ei — 
0). The lowest curve is one possible intermediate channel. 



higher partial wave contributions, while they exist, only 
contribute to rates at the fraction of a precent level. 

Because the electric field remains zero, parity is still 
a rigorously good quantum number. Therefore states 
of the same parity are not directly coupled. Nev- 
ertheless, the dominant loss channels in a magnetic 
field share the parity of the incident channel wave 
function, \i) = |22,-)|22,-)|00). Figure (a) illus- 
trates this by showing the total (black) and partial 
(color) inelastic rates as a function of the magnetic field. 
The loss rates shown correspond to the exit channels 
|/} = |10-)|22-)|22) (green), 1 11-) 1 22-) |41) (blue) , and 
|10— >|10— )|44) (red). 

Since direct coupling to these final channels is forbid- 
den to the dipolar interaction, all coupling must occur 
through some intermediate channel \int). Moreover, ow- 
ing to the parity selection rules in the matrix element 
l|15p. this intermediate channel must have parity quan- 
tum numbers e± = £2 = +• Since this coupling is second 
order, the dominant exit channels can consist of both 
d-wave (If — 2) and g-wave (If =4) contributions, in 
contrast to the electric field case. 

The primary feature of the inelastic rates in Fig. [S] 
(a) is that they decrease significantly at large field. This 
decline is the main reason for optimism regarding evap- 
orative cooling strategies in OH; an applied bias field of 




E-E f (K) 

FIG. 6: Inelastic rate constants for the 3-channel model sys- 
tem, Eq. 1171 as function of initial and final threshold separa- 
tion. The two curves are for d and (/-wave exit channels black 
and red, respectively. In g-wave channel if| evolves more 
slowly as the thresholds are separated. 

3000 Gauss can reduce the inelastic rate constant to be- 
low 2 x 10~ 13 cm 3 /sec, see Figure Q] The cause of this 
decrease can be traced directly to the relative separation 
of the incident and final channel thresholds, along with 
the indirect nature of the coupling. 

To see this, we reduce the model to its essential in- 
gredients: (1) a strong dipole-dipole interaction, (2) the 
relative motion of the thresholds as the magnetic field is 
tuned, (3) an extremely exothermic intermediate channel 
and (4) the centrifugal barrier in the final and interme- 
diate channels. The Hamiltonian for a reduced model is 
H Model = T + Vuodeh where T is the kinetic energy 
operator and VModei in matrix form is 

/ Ei a/R 3 \ 

V M odei=[ E f + Cf/R 2 b/R 3 .(17) 
\a/R 3 b/R 3 E lnt +c lnt /R 2 J 

Here Cj is a centrifugal repulsion Cj = h 2 lj(lj + l)/2m r , 
a and b are dipole-dipole coupling strengths, and Ej are 
the threshold energies for the j th channel. The channels 
{%, f, int} have quantum numbers (e^Y = (^1^2) = 

( ) and (£162)™* = (++)• The incident channel has 

partial wave i, = 0, while dipole coupling selection rules 
allow li n t = 2, and If = 2 or If = 4. 

The model Hamiltonian l|17fl explicitly excludes direct 
coupling between incident and final channels, whereas 
coupling is mediated through the int channel. Parame- 
ters characteristic of the physical problem are a = 0.12 
(a.u.), b = 0.10 (a.u.), E { =0,E f = -0.003 -> -0.1(if), 
and E mt = -0.17 (K), k = 0, If = 2 or l f = 4 and 
hnt — 2. Because of the energy gap separation losses 
to the intermediate channels are negligible. We find, in 
addition, that moving Ei nt has little effect on the rate 
constants for loss to channel /. 

Fig. shows the inelastic rates computed within this 
model. This three channel model does a reasonable job 
of mimicking the prominent features of the full calcula- 
tion, including the eventual and lasting decrease in rates 
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as the states are separated in energy. In addition, the g- 
wave rates decay more slowly as a function of field than 
do the d-wave rates, consistent with the full calculation 
(compare Fig. 0). The declining values of the rate con- 
stant cannot, however, be attributed to a simple overlap 
integral of the form Ijl6|l . since the incident and final 
channels are not directly coupled. We therefore present 
a more refined adiabatic analysis of this process in the 
next subsection. 



C. Adiabatic Analysis of the magnetic field case 

To understand the system's magnetic field behavior we 
analyze the reduced channel model (fTTjl in the adiabatic 
representation pHl2^|. This representation assumes that 
R is a "slow" coordinate. At every R we diagonalize the 
Hamiltonian in all remaining degrees of freedom. Since it 
is not rigorously true that R varies infinitely slowly, the 
residual nonadiabatic couplings can be accounted for in 
the kinetic energy operator. Written more formally we 
diagonalize 

K=(^ + v_ { n ) + Z z ) ( 18 > 

where the terms are the centrifugal barrier, potential ma- 
trix including dipole-dipole interaction and the Zeeman 
Hamiltonian. Diagonalizing the matrix in Eq. 1181 we 
get W\a(R)) = U a (R)\a(R)) where U a (R) are the eigen- 
values and \a(R)} are the eigenvectors. With the eigen- 
vectors we are able to form a linear transformation X(R) 
which transforms between the diabatic and adiabatic rep- 
resentations, ie X T W (R)X = U_(R). The eigenvalues 
and eigenvectors have radial dependence, but for nota- 
tional simplicity (R) will be suppressed hereafter. 

To distinguish between adiabatic and diabatic repre- 
sentations we use Greek letters (a, 13,...) to denote the 
adiabatic channels and Roman letters (i,j,...) to denote 
diabatic channels. When considering specific inelastic 
processes in the diabatic basis we denote initial and final 
channels as i and / and for the adiabatic channels as u 
and ip. In the limit R — » oo, the two sets of channels 
coincide. 

A partial set of adiabatic potential curves generated 
in this way is shown in Fig. [7] exhibiting an avoided 
crossing at R = 150. Thus molecules incident on the 
uppermost channel scatter primarily at large values of R. 
This point has been made in the past when an electric 
field is applied here we note that it is still true in 
zero electric field, and that scattering calculations can 
proceed without reference to short-range dynamics. 

The transformation between the representations is R 
dependent implying that the channel couplings shift from 
the potential to the kinetic energy operator. Using the 
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■0.02 -I , 1 , 1 , , , 1 , 1 

100 200 300 400 500 

R (a.u.) 




R(a.u.) 

FIG. 7: The relevant adiabatic potential curves for the OH- 
OH system. Shown are two different values of the final thresh- 
old energy E t — Ef = 5mK (a) and Ei — Ef = 15mK (b) 

adiabatic representation changes Eq. ^| to 

Here f = X T ^p. 

To get the channel couplings in the adiabatic picture 
we need matrix elements of the derivative operators, de- 
fined as P a /3 = (a\j^\P). We evaluate the the P a p ma- 
trix, the dominant off-diagonal channel coupling, using 
the Hellmann-Feynman theorem [22( 

p Efc; xl k vv kl x lf3 

Pa " ~ u a ~u (i ■ (20) 

Scattering amplitudes are then easily estimated in 
the adiabatic distorted-wave Born approximation (AD- 
WBA). Namely, we construct incident and final radial 
wave functions tf> ljV that propagate according to the adi- 
abatic potentials U L . V . In terms of these adiabatic wave- 
functions, the scattering T-matrix is given by an overlap 
integral analogous to Eqn. ^| 

nh 2 d/dRg+^/dR 
T llfi = {4> L \ -= 1^). (21) 
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FIG. 8: Inelastic rate constants as estimated by the adiabatic 
distorted wave Born approximation for the 3-channel system. 
The black curve is for a d-wave exit channel and the red for 
a g-wave exit channel. 



Here d/dR (d/dR) is the radial derivative operator act- 
ing to the left (right). The cross section for identical 



bosons is a Lip — 



87r |T t J 2 . 



From here we are able to 



numerically calculate a rate constant for inelastic loss 
j^mei _ Vl (x where v L is the asymptotic velocity given 
by y/2U L {R -+ oo)/m r . 

The result of the ADWBA is shown in figure El The 
two curves are for d (black) and g-wave (red) inelastic 
channels. Several key features are present that also oc- 
cur in the full calculation, namely: (1) the inelastic rate 
goes down with increasing threshold separation. (2) there 
is a zero in the rates as seen in Fig. (3) the g-wave in- 
elastic rate goes more slowly than the d-wave as seen in 
the model and the full calculation. The ADWBA ac- 
counts for all of these. The first feature, diminishing 
rates, still arises from an energy gap suppression, since 
the de Broglie wavelengths of incident and final channel 
still do not match well. In the ADWBA this process is 
further helped along by the fact that the residual channel 
coupling, represented by P, is localized near the avoided 
crossings of the adiabatic potential curves. 

The ADWBA helps to visualize this suppression, as 
shown by the sample wave functions in Fig. This fig- 
ure shows ip L , -j^ipif, and P Ltp for various values of Ef. 
Varying Et mimics the threshold motion of the system. 
The values of Ef of Figure are (a) Ef = —6mK, (b) 
E f = ~22mK, and (c) Ef = -62mK. The effect of 
the different EfS leave tpt, mostly unchanged, however tp^ 
becomes more exothermic and therefore more oscillatory 
(Xdb clearly shortens). Moreover, the dominant coupling 
region, where P Ltp peaks, moves to shorter R as Ef in- 
creases. This motion is obvious from the avoided crossing 
in Fig. □ 

The transition amplitude in the ADWBA is propor- 
tional to the integral of the product of the three quanti- 
ties in Fig. |21 Because of the shortening of the de Broglie 
wavelength in the exit channel, this integral will eventu- 
ally vanish, accounting for the zero in the inelastic rates. 
The total rate will, in general, not vanish, since there 
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200 300 400 500 

R(a.u.) 

FIG. 9: Illustrations of the origin of the zero in the partial 
rates. Each panel shows curves ip L , -r^tpip X 100 and P Llp xlO 9 . 



dR ^•fi ■ 

The plots are for different final energies and inelastic rate from 
the Born approximation are (a) Ef = —6mK and K' L ^ el = 
5 x 10~ n (b) E f = -22mK , K\ n / = 5 x 10" 14 (c) E f = 
-62mA', Kl™ el = 6 x 10" 13 , see text for details. 



are many exit channels, and they will experience the de- 
structive interference at different values of the threshold, 
hence at different fields. 

Finally, the g-wave inelastic rates are not so strongly 
affected by the separation of Ei and Ef because the g- 
wave centrifugal barrier is larger, meaning a greater en- 
ergy is required to change the wave function at short 
range such that velocity node can pass through the cou- 
pling region. The zero in this rate constant will thus 
occur at larger threshold separations. 
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V. CONCLUSIONS 

We have explored the influence of a magnetic field 
on the cold collision dynamics of polar molecules. The 
dipole-dipole interactions remain significant even in the 
absence of an electric field that polarizes the molecules. 
In general this implies that molecular orientations are un- 
stable in collisions, making magnetic trapping infeasible. 
We have found, however, that a suitably strong magnetic 
field can mitigate this instability. 

Beyond this result, we note that laboratory strength 
fields can exert comparable influence on cold collisions, if 
applied separately. A useful rule of thumb in this regard 



is that an electric field of 300 V/cm acting on a 1 Db 
dipole moment caues roughly the same energy shift as a 
100 gauss field acting on a 1 Bohr magneton magnetic 
moment. This raises the interesting question of how the 
two fields can be applied simultaneously, to exert even 
finer control over collision dynamics. This will be the 
subject of future investigations. 
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